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I. INTRODUCTION 

Spectroscopic [T] observations of absorption lines in 
the intergalactic medium (IGM) can be used to study 
the origin of large-scale structure in the universe, the 
history of star and galaxy formation, the metagalactic 
radiation field, and the chemical evolution of the IGM 
El |4] . Observations of lines from C IV, N V, O VI, and 
Si IV are routinely employed for these studies. These ob- 
servations provide important constraints for IGM studies 
when coupled with calculations of the ionization balance 
using codes such as CLOUDY [5]. However, the accu- 
racy with which one can infer the properties of the IGM 
is limited by uncertainties in the underlying atomic data. 

Of particular importance are reliable electron-ion re- 
combination data for the process known as dielectronic 
recombination (DR). This is the dominant recombina- 
tion process for most atomic ions under IGM conditions. 
Recently, Savin [2] has investigated the importance of 
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DR for CIV, NV, O VI, and Si IV. (Here, the convention 
of identifying each ion by its charge state before the re- 
combination process is used.) His work has shown that 
uncertainties in the DR data for these four ions limit our 
ability to constrain the metagalactic radiation field and 
the initial mass function for the earliest generations of 
stars. 

In the last several years, a series of measurements 
has been carried out to produce accurate DR data for 
CIV [6], NV [7 , OVI R, and sodium-like Si IV which 
will be presented in this work. Other experimental re- 
sults for DR rate coefficients of Na-like ions have been 
published for Si IV [9J, for FeXVI [10] and for Ni XVIII 
[TT] . Recently detailed calculations of the low energy DR 
resonance structure have been performed for Si IV [T2] 
within the framework of relativistic many-body pertur- 
bation theory (RMBPT). 

In the present work, experimental results for the Si IV 
recombination rate coefficient are presented. These 
were obtained employing the electron-ion merged-beams 
method at a heavy-ion storage ring. Orban et al. [S] 
measured the Si IV recombination rate coefficient in the 
electron-ion collision energy range 0-20 eV that com- 
prises DR resonances associated with 3s — > 3p and 
3s — > 3d (AN — 0) core excitations. Here an extended 
energy range of up to 186 eV was experimentally investi- 
gated. This additional range includes DR resonances as- 
sociated with 3s — > nl (n > 4), 2p — > nl (n = 3, . . . , oo) 
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(AN = 1,2,...) and 2s nl (n = 3,4) (AN = 1,2) 
core excitations. The corresponding excitation energies 
are listed in Table |T] except those for 2s — > nZ excitation. 
The present experimental Si IV merged-beams recombi- 
nation rate coefficient thus benchmarks theory for light, 
low charged sodium-like ions over a wide range of ener- 
gies. 

TABLE I: Excitation energies -E exc for the excitation of Si IV 
2s 2 2p 6 3s 2 Si/2 ground state to 2p 6 nl and to 2p 5 3s nl (both 
n — 3,4:) states that are relevant in the present work. For the 
derivation of the 2p 5 3s nl energies, results from experimental 
Auger spectroscopy [13] were added to the Si IV ionization 
energy of 45.14179 eV [14]. 



excited state -Eexc (eV) Reference 

2p 6 3p 2 Pi /2 8.839 [T5] 

2p 6 3p 2 P 3/2 8.896 [15] 

2p 6 3d 2 D 6/2 19.884 [14] 

2p 6 3d 2 D 3/2 19.884 [14] 

2p 6 4s 2 Si /2 24.050 [H] 

2p 6 4p 2 P 1/2 27.062 [14] 

2p 6 4p 2 P 3/2 27.082 [14] 

2p 6 4d 2 D 5/2 30.997 [14] 

2p 6 4d 2 D 3/2 30.997 [H] 

2p 6 4/ 2 F 5/2 31.508 [14] 

2p 6 4/ 2 F 7/2 31.508 QI] 

2p 5 3s 2 2 P 3/2 99.06 [T3] 

2p 5 3s 2 2 P 1/2 99.68 [T3] 

2p 5 3s3p( 3 P) 4 S 3/2 104.8 [13] 

2p 5 3s 3p ( 3 P) 4 D 106.2 [13] 

2p 5 3s 3p ( 3 P) 4 P 106.9 [13] 

2p 5 3s3p( 3 P) 2 D 107.5 [T3] 

2p 5 3s3p( 3 P) 2 Si /2 108.4 [13] 

2p 5 3s3p( 1 P) 2 D 5/2 110.9 [T3] 

2p 5 3s3p( 1 P) 2 D 3/2 111.5 [13] 

2p 5 3s3p( 1 P) 2 P 112.0 [13] 

2p 5 3s3p( 1 P) 2 Si /2 112.7 [13] 

2p s 3s3d( 3 D) 4 P 3/2 119.2 [13] 

2p 5 3s 3d ( 3 D) 4 F 119.9 [13] 

2p 5 3s3d( 3 D) 4 D 7/2 121.1 [13] 

2p 5 3s3d( 1 D) 2 F 7/2 124.1 [13] 

2p 5 3s3d( 1 D) 2 D 3/2 124.6 [13] 

2p 5 3s3d( 1 D) 2 D 1/2 125.1 [13] 

2p s 3s( 3 P)4s 4 P 3/2 125.8 [T3] 

2p 5 3s( 3 P)4s 2 P 126.3 [13] 

2p 5 3s( 1 P)4s 2 P 3/2 126.9 [T3] 

2p 5 3s4p( 3 P) 4 S 128.4 [13] 

2p 5 3s4p( 3 P) 2 P 129.5 [13] 

2p 5 3s4p( 1 P) 2 P 131.0 [13] 

2p 5 3s4d( 3 D) 4 F 133.3 [T3] 

2p 5 3s4d( 3 D) 2 P 134.7 [T3] 



In this work, the experimental results are compared 
to theoretical results using the multiconfiguration Dirac- 
Fock (MCDF) method, in particular at low electron-ion 
collision energies. At these energies, the calculation of 
accurate resonance positions and strengths is extremely 
critical for the derivation of reliable plasma DR rate coef- 
ficients. Therefore, the present work is partly an investi- 



gation into the capabilities of the MCDF method for cal- 
culating accurate DR resonance energies and strengths. 
The MCDF method is an ab-initio method that is ap- 
plicable to target atoms and ions with an arbitrary shell 
structure. In contrast, the RMBPT is currently limited 
to quasi one-electron target ions, owing to its inherent 
complexity. For such systems, however, it usually yields 
more accurate results than any other theoretical method. 
Below, a detailed comparison between the present MCDF 
results and the RMPBT results of Orban et al. [12] for 
DR of sodiumlike Si IV will be presented. 

The present paper is organized as follows. The theo- 
retical procedure is outlined in Section [TTJ. In Section [III] 
the experimental procedure is described. Experimental 
and theoretical results are presented and compared in 
Section |IV[ The Si IV to Si III experimental and theo- 
retical DR rate coefficients in a Maxwellian plasma are 
derived and compared to recent theoretical and experi- 
mental results in Section[V] Conclusions will be presented 
in Section [VI] 

II. THEORY AND IMPLEMENTATION 

In most computations, DR of an A/"-electron target ion 
in the initial state \i) is handled as a two-step process 
in which first an electron is captured resonantly from the 
continuum forming an (J\f + l)-clcctron state \d) with 
the captured electron and one of the target core electrons 
both now in an excited level. In a second step, this doubly 
excited state cither decays radiatively by the emission of 
one or more photons to some final state |/) which lies 
below the ionization threshold of the ion, or it returns 
by Auger electron emission to the initial ionization stage 
of the target. If interference between the radiative and 
nonradiative capture of the electron in the field of the 
target is negligible, the integrated DR cross section for 
an isolated resonance, the so-called resonance strength, 
can be expressed in terms of the Auger and radiative 
rates of the intermediate state \d) by 

/oo 
<j dr (E) dE 
-OO 

_ 2ir 2 h A a (i -> d) A r (d -> /) 

where ki denotes the wave number and E the energy of 
the incident electron, Yd the width of the doubly-excited 
state, and A a (i — > d) is the rate for the dielectronic cap- 
ture from the initial into the doubly excited state. Using 
the principle of detailed balance, the dielectronic cap- 
ture rate A a (i — > d) is equal to [gd/di A a (d — > i)] where 
A a (d — * i) is the Auger rate, and gi and gd are the sta- 
tistical factors of the initial and intermediate states, re- 
spectively. A r (d — * /) refers to the rate for the radiative 
stabilization to the state |/). The strength as defined by 
Equation Q refers to the area under the DR cross section 
curve. Hence, it is usually given in units of cm 2 eV. The 
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total width Yd is determined by all possible decay chan- 
nels of the resonant state \d) and is given, in first-order 
perturbation theory, by 



(2) 



where the sums are over all the individual Auger and ra- 
diative rates (widths) of the intermediate state \d). The 
use of the resonance strength S(i — > d — » /) is appropri- 
ate for resonances whose width Td is small compared to 
its energy position E Tes = Ed — Ei where Ed and Ei are 
the energy of the doubly excited state and initial state of 
the target ion, respectively, and if the energy dependent 
DR cross section has a Lorentzian profile [TFj 



T d /2 



7T {E ICS -EY + T 2 JA 



(3) 



Owing to the energy of the incident electron [cf . the l/kf 
factor in S\, the resonance strength increases rapidly to- 
wards the DR threshold and makes the process sensitive 
to low kinetic energies of the incoming electrons. 

At storage rings, the radiative stabilization is often 
not observed explicitly and, hence in Equation ([!]), the 
radiative rate for the individual transition d — > / has to 
be replaced with 



A r (d^f) 



f" 



/")> 



where the summation extends over all states \f") below 
the ionization threshold that can be reached by radia- 
tive transition from the doubly excited state \d). As 
seen from Equation ([I]), the observed DR strength for 
a given resonance is proportional to the capture rate 
A a (i — * d) and the total radiative rate. For most light 
and medium elements, moreover, the magnitude of the 
resonance strengths is often determined by the radiative 
decay since the emission of photons from the doubly ex- 
cited state is then much slower than its autoionization 
and, hence, A a A r /(A a + A r ) ss A r . 

The doubly excited state \d) is often one out of a large 
number of highly correlated states embedded in the con- 
tinuum of the target ion. For these resonances, special 
care has to be taken in calculating both the individual 
as well as total rates accurately. To describe the ground 
and excited state of multiply charged ions, the MCDF 
method has been found to be a versatile tool for the com- 
putation of the many-electron energies and decay prop- 
erties, especially if inner-shell electrons or several open 
shells are involved in the computations [TTj HH] . In the 
MCDF method, an atomic state is approximated by a 
linear combination of so-called configuration state func- 
tions (CSF) of the same symmetry 



MPJ) = £ Cr{a) \ lr PJ) , 



(4) 



where n c is the number of CSF, {c r (a)} denotes the rep- 
resentation of the atomic state in this CSF basis, P is the 
parity, J is the total angular momentum, and j r is a set of 
quantum numbers for a unique specifying of the many- 
electron basis states. In most standard computations, 
the CSF are constructed as antisymmetrized products of 
a common set of orthonormal orbitals and are optimized 
on the basis of the Dirac-Coulomb Hamiltonian. Further 
relativistic contributions to the representation {c r (a)} of 
the atomic states are then added, owing to the given re- 
quirements, by diagonalizing the Dirac-Coulomb-Breit 
Hamiltonian matrix in first-order perturbation theory. 
For multiply charged ions, an estimate of the dominant 
QED contributions (i.e., the self-energy and vacuum po- 
larization of the electronic cloud) might be taken into 
account by means of scaled hydrogenic values. But QED 
plays a negligible role for Si 3+ ions in the present analy- 
sis. 

Obviously, the calculation of the radiative and nonra- 
diative decay rates of the intermediate resonances \d) is 
central to the computation of any DR spectrum which is 
to be compared with experiment. In the MCDF model, 
both the radiative as well as the Auger matrix elements 
are derived from computation of the corresponding inter- 
action matrix within the CSF basis using the expansion 
@. To deter mine the rates, the wave functions from 
the GRASP 92 [TH] and RATIP gU] codes have been ap- 
plied in the present work which allows one to incorporate 
both the dominant relativistic and correlation effects on 
the same footings. However, since the computation of 
the transition probabilities has been considered at many 
places elsewhere [THUH], we shall mention here only that 
the radiative rates 



A r (d^f) = 



3 c 2 2 J d + 1 

Y,\^f{JfPf)\H- t ^L)\MJdP d )) 
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where w^j is the emitted photon's frequency, are asso- 
ciated with the reduced matrix elements of the interac- 
tion H 7 (ttL) of the atomic electrons with the multipole 
components of the radiation field |21j , where it = and 
7r = 1 refer to electric and magnetic multipoles, respec- 
tively, with the multipolarity L. For light and medium 
elements, it is of course sufficient to include the electric- 
dipolc (L = 1, 7r = 1) decay while the contributions 
from higher multipoles to the radiative stabilization re- 
main negligible. To obtain the total radiative rate, i.e., to 
include the summation over all lower states |/') in Equa- 
tion ([2]) , often a large number of transition rates have to 
be compiled with similar accuracy. For this reason the 
computation of the DR spectrum at higher energies may 
become difficult. In the RATIP program, the transition 
probabilities are computed by the REOS component [25] 
including, if appropriate, the rearrangement of the elec- 
tron density in course of the decay. The explicit consid- 
eration of the electron density's rearrangement was not 
included here, because this would require separate calcu- 
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lations for the intermediate resonance states and the final 
states of the radiative stabilization. In the present calcu- 
lations, however, both sets of states were always treated 
together in order to keep the computations feasible. 

Unlike the treatment of radiative stabilization, calcu- 
lation of the Auger rates requires the coupling of the 
bound-state electrons (of the target ion) to the electron 
continuum. If, for the sake of simplicity, we neglect the 
interaction between different Auger channels (i.e. within 
the continuum), the Auger rates are given by 

A a {d -» j) = 2tt x (6) 
2 | (WjiJjPj), en c ) J d P d \H-E d \ MJdPd)) | 2 , 

where H is the atomic Hamiltonian, E d the total energy 
of the doubly excited state and e — E d — Ejthe kinetic 
energy of the emitted electron. In Equation pi, the sum- 
mation over k c = ±(j c +l/2) for I = j c ±l/2 , where I and 
j c are the angular momentum and total angular momen- 
tum, respectively, of the incoming or outgoing electron, 
extends over all partial waves of the outgoing electron 
which can be coupled to the target ion state \ipj(Jj,Pj)) 
with the condition to conserve the total angular momen- 
tum J d and parity P d of the intermediate state. If, more- 
over, a common set of orthonormal orbitals is used for 
the representation of the intermediate state \d) and the 
final ionic states then the operator (H — E) « V 
can be replaced by the electron-electron interaction op- 
erator. For most light and medium elements, it is again 
sufficient to include the instantaneous Coulomb repulsion 
between the electrons but to omit the relativistic Breit 
contributions as they were found important only for the 
Auger emission of highly-charged ions [23 , 24 , 25] [26] . 

The restriction of the electron-electron interaction in 
the computation of the matrix elements in Equation ([6]) 
is common practice, even though the orbital functions 
for the doubly excited state \d) and the final ionic state 
|j) are not quite orthogonal to each other. This treat- 
ment has been implemented therefore also in the AUGER 
component of the RATIP program, in which the con- 
tinuum spinors are solved within a spherical but level 
dependent potential of the final ion (the so-called optimal 
level scheme in the GRASP92 program). This scheme 
also includes the exchange interaction of the emitted elec- 
tron with the bound-state density. Often, the number of 
the possible scattering states \(ipj(JjPj), en c ) J d P d ) of a 
system increases rapidly with the number of intermediate 
ion states as the free electrons may couple in quite differ- 
ent ways to the bound-state electrons. For further details 
on the computations of the Auger matrix elements and 
relative intensities, we refer the reader to [TB] and [27]. 

Apart from the individual rates, any helpful calcula- 
tion of DR spectra critically depends on the proper con- 
trol and handling of the various decay branches which, at 
least in principle, should be treated on the same basis in 
order to obtain a consistent spectrum. This need for an 
efficient handling of the decay channels concerns those 



open-shell configurations in which one-particle excita- 
tions with large principal quantum numbers 'mix into' 
the low-energy part of the DR spectrum. Such a mix- 
ing occurs especially for the excitation of electrons be- 
neath the valence shell whose energies are comparable 
with the valence excitation into Rydberg states. There- 
fore, in order to facilitate this handling of the various 
decay branches, a new component (DIEREC) has been 
developed recently in the framework of the RATIP pro- 
gram. This component now supports both, the compu- 
tation of individual S(i — > d — > f) and total resonance 
strengths S(i — > d) and enables us to simulate the low- 
energy DR spectra as observed by experiment. Inter- 
nally, of course, this new component makes use of the 
REOS and AUGER components from above and allows 
us, if necessary, to incorporate also higher multipoles 
in the computation of the (total) radiative rates. Be- 
cause of the finite-difference method of the GRASP92 
and RATIP programs and owing to the associated size 
of the configuration expansions, however, excitation to 
states with principal quantum numbers n > 8 can of- 
ten not be treated by the code explicitly. Atomic struc- 
ture codes, such as, e.g., AUTOSTRUCTURE [2B], that 
are geared towards the generation of DR rate coefficients 
for plasma physical applications, treat high-n electrons 
within a hydrogenic approximation. This has not yet 
been implemented in the RATIP program suite. 



III. EXPERIMENT 

The experiment was performed at the heavy-ion stor- 
age ring TSR of the Max-Planck-Institut fur Kernphysik 
(MPI-K) in Heidelberg, Germany. Details of the merged- 
beams technique using the MPI-K electron cooler have 
been described previously [6] |2"9l |3tJl I3T] I32] . 

In the present experiment a beam of 28 Si 3+ was pro- 
vided by the tandem accelerator of the MPI-K linear ac- 
celerator facility at an energy of about 1.1 MeV/u. The 
ion beam was injected into the storage ring where it was 
collinearly overlapped with the cooler electron beam. For 
efficient electron cooling, the electron velocity has to be 
close to the ion velocity in the ring. To fulfill this velocity- 
matching condition, the laboratory energy of the elec- 
trons was set to the 'cooling' energy E c ss 593 eV. The 
beam current was accumulated by multi-turn injection 
and 'ecool stacking' [35]. Ion currents of typically 10- 
50 fiA were reached. 

In contrast to previous experiments, where the electron 
beam of the cooler was also used as an electron target 
for recombination experiments, in the present experiment 
a newly installed separate electron beam [33] was used. 
This additional electron beam is hereafter denoted as the 
electron target. As in the electron cooler, the electron 
beam of the electron target is also guided by a magnetic 
field and overlaps the ion beam in a straight section of 
£3 1.5 m length. 

Conceptually, the experimental procedures for measur- 
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ing recombination rate coefficients with the electron tar- 
get are the same as those applied previously with the 
electron cooler. However, there are advantages when a 
separate electron target is employed for recombination 
measurements. First, the electron cooler can be used con- 
tinuously for the cooling of the ion beam. Thus, the low 
velocity and spatial spread of the ion beam is maintained 
at all times. Second, the electron target was specifically 
designed for providing an electron beam with a very low 
initial energy spread [34] . Both advantages yield a higher 
experimental resolution in the present measurement as 
compared to previous measurements using only the elec- 
tron cooler. 

After injection into the storage ring, the ions were 
cooled for a few seconds before the recombination mea- 
surements started. Recombined ions were separated from 
the circulating beam in the first dipole magnet down- 
stream of the electron target and counted by a sin- 
gle particle scintillation detector with nearly 100% ef- 
ficiency During the recombination measurement, the 
electron energy of the electron target beam was alternat- 
ingly chopped between measurement (E m ) and reference 
(E r ) energies by switching the acceleration voltage for 
the target electron beam accordingly. At the same time 
the electron cooler was held constant at cooling energy. 
The reference measurement was made to determine de- 
tector background. Therefore, the reference energy E r 
was chosen in a range of the spectrum where no DR res- 
onances occurred. The measurement and reference in- 
terval of data collection were 10 ms each. Between the 
voltage jump and the data collection interval there was a 
settling time of 5 ms to allow the power supplies to reach 
the preset values. The data were collected in overlapping 
data sets ranging from laboratory energies E m = 561 eV 
to E m = 1420 eV. The chopping pattern was E!^-E r - 
E^-Er with v = 1,2,3.... With each step v in the 
chopping pattern E m was changed by 0.15 eV in the lab- 
oratory frame whereas E r was kept fixed. The merged- 
beams rate coefficient is derived from the background 
subtracted recombination count rate using [29] [35] 



[R(E m ) — R(E m , E r )\ 
(1 - piPJrpiciEML/C 



a{E r ) 



n c {E r ) 
n e (E„ 



(7) 

where R is the recombination count rate, rj is the de- 
tection efficiency, n c is the electron density in the in- 
teraction region, Ni is the number of ions in the ring 
and C — 55.4 m is the circumference of the TSR. In 
a first approach the nominal length L of the electron- 
ion overlap region in the electron target section is set to 
L — 1 .476 m which is the length of the solenoid providing 
the axial magnetic field along the interaction region. For 
the detailed analysis we applied a toroid correction [30] 
that accounts for the contribution of the merging and 
de-merging sections in the toroidal magnetic fields of the 
electron target. The ion velocity and the electron veloc- 
ity are Vi — (3\c and v c — (3 c c, respectively, where c is the 
speed of light. Under the present conditions, the factor 
1/(1 — A/?e) deviates by less than 0.4% from unity. The 



second term in Equation ^ is a small correction that 
re-adds the recombination rate coefficient a(E r ) at the 
reference point. This contribution was calculated using a 
hydrogenic formula [5] for nonresonant radiative recom- 
bination (RR) . The systematic experimental uncertainty 
of the merged-beams rate coefficient is estimated to be 
±18% at lex confidence level. This uncertainty stems 
mostly from the ion current measurement (±15% includ- 
ing also the error of the toroid correction) [SD] and the 
determination of the electron density (±10%). 

The experimental electron energy distribution is 
best described as a flattened Maxwellian distribution 
which is characterized by the longitudinal and trans- 
verse temperatures Tji and T±. The experimental 
energy resolution is approximately given by A£ = 
[(ln(2) /c B Ii) 2 + 161n(2) EksT^} 1 / 2 , where E denotes the 
relative electron-ion energy. For the reduction of Tj_ the 
target electron beam was adiabatically expanded |31j by 
factors of up to 28. 

For the measurements we used two different cathodes, 
a thermionic cathode and a photocathode [3B] . The elec- 
tron current produced by the thermionic cathode was 
4 mA and that of the photocathode 0.5 mA. 

The cryogenic photocathode provides electrons with a 
laboratory energy spread of about 10 meV [35]. The ex- 
pansion factor was 20. The electron beam temperatures 
were estimated by fitting a simulated (Sec. IV A I merged- 
beams recombination rate-coefficient comprising contri- 
butions by DR and RR to the measured spectrum. A 
transverse temperature of k^T± « 0.9 meV and a longi- 
tudinal temperature of fcaTy ~ 35 /j,eV were found. With 
the above given temperature values, the experimental en- 
ergy spread in the center of mass frame thus amounts to 
AE = 6 meV at E = 0.1 eV, AE = 20 meV at E = 1 eV, 
AE = 62 meV at E = 10 eV, and AE = 197 meV at 
E = 100 eV. 

The temperature of the thermionic cathode is typi- 
cally about 1300 K. The expansion factor was 28. The 
temperatures were again derived by fitting a simulated 
spectrum to the measured one. A transverse tempera- 
ture of fc B ^J_ ~ 4.4 meV and a longitudinal temperature 
of /caTj| ~ 100 /xeV were determined. With these tem- 
peratures, the experimental energy spread thus amounts 
to AE = 11 meV at E = 0.1 eV, AE = 33 meV 
at E = 1 eV, AE = 105 meV at E = 10 eV, and 
AE = 333 meV at E = 100 eV. 



IV. MERGED-BEAMS RATE COEFFICIENTS 

A. Experimental results 

Figure [T] shows the measured merged-beams recombi- 
nation rate coefficient in the energy range 0-9.2 eV. For 
this measurement the photocathode was used. In this 
range all DR resonances associated with 3s — > 3p core 
excitations appear. The Rydberg series of 2p 6 3pn'l' DR 
resonances converges to its series limit at Eoo = 8.877 eV 
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FIG. 1: Measured Si IV to Si III merged- beams electron-ion 
recombination rate coefficient in the energy range dominated 
by DR resonances associated with 3s — > 3p core excitations. 
The vertical bars below the spectrum denote the 2p 6 3p n'l' 
DR resonance positions as expected on the basis of the hydro- 
genic Rydberg formula (Eq. Rjl. Note that resonances up to 
n = 15 are individually resolved. The vertical arrows around 
3 eV denote the positions of the 2p 6 3d 2 3 F and 2p 6 3d 2 1 G 
doubly excited states. The inset shows the recombination 
rate coefficient at energies below 2 eV, where it is dominated 
by 2p 6 3p 41' resonances. 



|15j . This value is the weighted average of the 3pi/2 
and 3p 3 /2 limits. The 3p fine structure splitting of 
0.057169 eV is not resolved. Additionally, the 2p 6 3d3l 
and 2p e 3d4s resonances associated with the 3s — > 3c? 
excitation are expected to appear at these low energies. 
As pointed out by Orban et al. jl2j, in this range the 
most notable 3s — > 3d contributions arc due to 2p 6 3d 2 3 F 
and 2p e 3d 2 1 G doubly excited states whose energy posi- 
tions are determined in our experimental data to 2.97 
and 3.13 eV, respectively. 

The experimental energy scale was fine-tuned by mul- 
tiplying the nominal electron-ion collision energies with 
an energy independent factor, so that so that the posi- 
tions of the 2p e 3pn'l' (7 < n < 10) resonances matched 
their calculated positions. For the high-n values the po- 
sition of the Rydberg resonances can be estimated from 
the hydrogenic Rydberg formula 



E„ — E„ 



R 



g 



(8) 



with the Rydberg constant R and the charge state q = 3 
of the initial Si IV ion. The calibration factor for the 
energy axis deviated from unity by less than 1%. 

The much less intense 2p 6 3d n'l' resonances with n' > 
4 associated with 3s — ► 3d core excitations are expected 
to occur in the energy range 9-20 eV. The recombina- 
tion rate coefficient measured in this energy range, using 
the thermionic cathode, is shown in Figure [2j Individual 
resonances of the associated Rydberg series are barely 
visible with rather large statistical uncertainties. Nev- 
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FIG. 2: Measured Si IV to Si III merged-beams electron-ion 
recombination rate coefficient in the energy range of DR reso- 
nances associated with 3s — ► 3d core excitations. The vertical 
bars denote the 2p 6 Zdn'l' DR resonance positions as expected 
on the basis of the hydrogenic Rydberg formula [Eq. pj]. The 
error bars on each data point show the statistical uncertain- 
ties. 



ertheless the 2p e 3dn'l' series limit at 19.884 eV [15] is 
clearly discernable. The maximum of the DR rate coef- 
ficient at the 2p 6 3dn'l' series limit is about two orders 
of magnitude smaller than that at the 2p 6 3pn'l' series 
limit. The energy range 20-69 eV is not shown because 
there is no significant structure exceeding 10~ 12 cm 3 s^ 1 , 
but in principle it comprises 2p 6 nl n'l' resonances asso- 
ciated with 3s — > nl excitation with n > 4. 




75 80 85 90 95 

Electron-ion collision energy (eV) 



FIG. 3: Measured Si IV to Si III merged-beams recombination 
rate coefficient in the electron-ion collision-energy region of 
resonances associated with AiV > 1 2s — » nl and 2p — > nl 
(n > 3) inner shell core excitations. 

In Figure [3] we show the measured merged-beams re- 
combination rate coefficient in the energy range 69- 
130 eV. This was measured using the thermionic cath- 
ode. The resonances in this energy range are associated 
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with A TV > 1 DR via 2p — > nl (n > 3) inner shell core ex- 
citations. Due to the large number of excitations in the 
displayed energy range (cf. Table [I]) it is prohibitive to 
assign individual resonances. The main contribution to 
the rate coefficient is most probably due to 2p — ► 3/ core 
excitations. Theoretical calculations suggest that DR by 
2s excitations is insignificant |37j . The energy range from 
130 up to 186 eV was also scanned and found to fluctuate 



with peak-to-peak variations up to ±10 12 cm 3 s 



not 



showing any significant structure. 

The calculated RR rate coefficient as well as 25 DR 
resonances fitted to to the measured DR spectrum in the 
energy region below 1.5 eV is shown in Figure [4] This 
fit is independent of the theoretical predictions and is ex- 



plained in detail in Section V A For the comparison with 
the theoretical calculation shown in the Figures [5j [6] and 
|7] the non-resonant RR contribution was subtracted from 
the measured merged-beams recombination rate coeffi- 
cient at all energies. The merged-beams RR rate coeffi- 
cient was derived by convolving the RR cross section with 
the experimental electron energy distribution. The RR 
cross section was calculated with a hydrogenic formula 
[Eq. (12) of Ref. [6], taking into account field ionization 
of loosely bound high Rydberg electrons inside the stor- 
age ring bending magnets Schippers et al. 
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FIG. 4: Measured Si IV to Si III merged-beams rate coefficient 
at energies below 1.5 eV (circles and grey shaded area). The 
solid curve is the sum of a fit comprising 25 fitted DR res- 
onances (dashed curves cf. Sec. VAl and the non-resonant 
rate coefficient due to RR (dotted curve). The inset shows 
the same curves, but in more detail in the energy range 1— 
100 meV. 



B. Theoretical results 

To simulate the observed DR spectra at different ener- 
gies of the initially free electron, detailed computations 
have been carried out for the 3s — > 3pn'V (n' = 3, . . . ,6) 
and 3s — > 3d n'V (n' = 3, ... ,6) resonant electron cap- 



ture and its subsequent radiative stabilization. The first 
DR resonance (3p4d 1 i?2) is found at about 0.1 eV, i.e., 
very close to zero energy of the incident electron. The 
exact energies and strengths of resonances close to the 
threshold are particularly important for the determina- 
tion of reliable plasma recombination rate coefficients at 
low electron temperatures. Therefore, special care was 
taken with regard to the doubly-excited magnesium-like 
states from the 31 n'V (n' > 3) final configurations of the 
recombined ion. 

For the 0-6 eV low-energy part of the DR spectrum, 
a series of computations has been carried out. In the 
first approach, we included all the levels within the 3/ 3V 
and 31 AV configurations of the recombined ion. Apart 
from these low-lying levels with energies both below and 
above the threshold, we incorporated in a further step 
also the levels of the 31 51' and 3/ 61' configurations as 
well as later all the nln'V levels with n < 5 and n' < 6, 
respectively. To keep the number of CSF manageable, 
levels with higher n and n! were not treated. As known 
from previous computations [38, 39 for the low-lying lev- 
els of multiply charged ions, such a systematic enlarge- 
ment of the wave function expansion Q typically im- 
proves the positions of the resonances significantly. For 
magnesium-like ions, moreover, many of the 31 n'V config- 
urations 'overlap' with each other in energy and, hence, 
'new' resonances may appear in the calculated low-energy 
part of the theoretical DR spectrum, if the configuration 
space is increased. Using the single and double excita- 
tions from above, we obtained an expansion of up to 1073 
CSF for the intermediate and the final-state wave func- 
tions of the Si III ions. 

For the 2p — > 3/ n'V part of the DR spectrum we could 
include only levels of 2p 5 3s 313V and 2p 5 3s3lAV , since 
the number of open shells involved is increased as com- 
pared to the 3s — > 31 n'V excitations. In the calculation 
the associated resonances appeared in the energy region 
69-94 eV of the incident electrons. For these inner-shell 
excited spectra, further contributions from the core po- 
larization or core-core excitations need to be omitted ow- 
ing to the size requirements of the corresponding wave 
function expansions. The incorporation of double excita- 
tions from the 2s and 2p shells would result in expansions 
of several hundred thousand CSF, i.e., a size which is 
unfeasible for the computation of DR and autoionization 
properties. 

Figure [5] displays the experimental Si IV DR spectrum 
in the energy region of 0-1.5 eV, with the nonresonant 
'background' subtracted from the experimental results. 
In this figure, the observed spectrum is compared with 
theoretical results from different approximations. In all 
of these calculations, our 'best' wave function expansion 
has been applied, including the single and double ex- 
citations as discussed above. The computations differ 
however in the set of the one-electron orbital functions 
used for the representation of the initial sodium-like ions. 
Figure |5ja), for instance, shows the spectrum in which 
both the initial and final states of the recombined ion 
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were described by a common set of orbitals, neglecting 
the rearrangement of the electron density in the course 
of the dielectronic capture (or decay) of the ions. Appar- 
ently, quite a strong effect arises from this rearrangement 
of the electron density as seen from Figure |5jb) and^c), 
for which two independent sets of orbital functions were 
utilized in the representation of the initial and the recom- 
bined ion states. Figures [5jb) and[5jc) only differ in the 
treatment of the exchange interaction for the incoming 
electron. While Figure |5jb) shows the simulation for a 
static potential due to the charge distribution of the ini- 
tial ion ('no exchange'), J5^c) incorporates the exchange 
interaction of the incident electron with regard to the 
bound-state density. Therefore, Figure [5jc) represents 
our best approximation within the MCDF approach. De- 
spite the fact that the orbitals are not quite orthogonal 
in the computation of the two-particle matrix elements, 
the Auger amplitudes [Eq. ([6])] were evaluated by using 
the techniques of Racah's algebra, i.e., for assuming or- 
thogonality for all inactive electrons in these transition 
amplitudes [40] . 

Figure |5^d) displays the comparison of the present ex- 
perimental data with the RMBPT result of Orban ct 
al. |12j convolved with the electron energy distribution 
of the TSR photocathode electron beam. Compared 
to their experiment, where electron beam temperatures 
k-eTn = 0.25 meV and k B T± = 10 meV were found, the 
energy resolution is higher in the present experiment. 
(fc B Tj| = 0.035 meV, k B T± = 0.9 meV, see Section [Til I . 
Our high-resolution experiment provides a more stringent 
test of the RMPBT calculation which represents the mea- 
surement almost perfectly but cannot easily be extended 
to energies beyond 1.4 eV where an increasing number of 
resonances can contribute to the DR spectrum. 

Table [n] displays the assignment and position of the 
28 lowest resonances in the energy region E ICS < 1.5 eV 
calculated with the MCDF method as well as the corre- 
sponding RMBPT results of Orban et al. [12]. In addi- 
tion to the energies of these resonances, this table also 
includes the weights of the dominant LS terms as well 
as the resonance strengths. The weights of the leading 
LS terms have been obtained by a unitary transforma- 
tion of the wave functions from the j j-coupled into a 
LS'-coupled basis |41j . While the lowest 11 resonances 
appear rather pure in LS coupling (> 92 %), some larger 
admixtures are found for a few of the higher-lying reso- 
nances. As mentioned above, all data in this table cor- 
respond to our best representation of the resonances and 
by including the effects of the rearrangement of the elec- 
tron density and the exchange interaction (cf. Fig. 
The lowest resonances in the 3s — > 3pn'l' part of the DR 
spectrum belong to the 3p4d 1 D% level, followed by the 
fine-structure levels of the 3p4o? 3 F° term. Apparently, 
all levels from the 3p 2 and 3p 3d configurations are below 
the DR threshold. 

At energies below 0.9 eV the RMBPT results of Orban 
et al. |12j are in excellent agreement with the present 
measurements [Fig. [5jd)]. At higher energies, where 
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Electron-ion collision energy (eV) 

FIG. 5: Comparison of the experimental Si IV DR spectrum 
(open circles) in the energy region — 1.5 eV with various the- 
oretical results (solid curves). Figures (a), (b), and (c) show 
our multiconfiguration Dirac-Fock (MCDF) method results 
and (d) published relativistic many-body perturbation theory 
(RMBPT) calculations (Ref. [12]). The nonresonant part of 
the recombination rate coefficient due to radiative recombina- 
tion (RR) was subtracted from the measured curve (see text). 
In the MCDF calculation three different approximations have 
been applied to the enlarged CSF basis including all the fine- 
structure states from the nl n'l' (n, n' = 3, 4, 5, 6) configura- 
tions: (a) using the orbital functions from the doubly-excited 
31 n'l' 2S+1 Lj levels for the representation of both, the ini- 
tial and final states; (b) with an independent optimization of 
the sodium- and magnesium-like states but without including 
the exchange interaction between the initially free electron 
and the bound-state electrons; (c) the same as in (b) but 
by incorporating the exchange interaction with regard to the 
bound-state density. 



the MCDF results reproduce the experimental findings 
slightly better than at lower energies, the RMBPT res- 
onance positions are at somewhat too high energies as 
can also be seen in the comparison with the experimen- 
tal data of Orban et al. [T2J. In our work, this slight 
discrepancy is more pronounced because of the increased 
experimental resolution and because of reduced statis- 
tical uncertainties in the present experiment. With a 
few exceptions, MCDF and RMBPT resonance energies 
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TABLE II: Comparison of the present multiconfiguration 
Dirac-Fock (MCDF) resonance parameters with the results 
of the relativistic many-body perturbation theory (RMBPT, 
Ref . [12] ) for all 3p \l and 3p 5s DR resonances associated with 
3s — > 3p core excitations. The listed quantities are resonance 
energies E [cs and strengths S [Eq. Q]. The dominant LS- 
terms of the resonance states are given in the first column 
along with their weights in the representation of the wave 
functions. The weights are derived from the MCDF calcula- 
tion. The states are listed in the order of increasing MCDF 
resonance energies. RMBPT resonance energies which appear 
in a different order [T^] are marked by an asterisk in the third 
column. 



dominating 


E TeB (eV) 


S (1CT 20 eV cm 2 ) 


LS term 


MCDF RMBPT 


MCDF RMBPT 



6.88 
3.77 
4.95 
5.58 
1.78 
2.96 
4.01 
1.42 
0.83 
0.27 
1.72 
3.20 
2.50 
1.46 
2.96 
0.24 
0.71 
2.48 
1.14 
3.13 
3.75 
2.23 
0.78 
1.01 
1.56 
1.13 
0.65 
0.54 



agree with one another to within ~ 50 meV, often even 
to within ~ 20 meV. 

Figures [6] and [7J display our theoretical MCDF DR 
spectrum compared with experiment for AN = and 
AiV = 1 DR. As described above, all nln'V configura- 
tions with 3 < n, nf < 6 have been taken into account 
for AiV = DR. The incorporation of further configu- 
rations with even higher principal quantum numbers n' 
has no effect upon the low-lying resonances for energies 
E rcs < 6 eV above the threshold. Although the basic fea- 
tures are well described in these spectra, some deviations 
in the positions and strengths of the peaks remain which 
we attribute to neglected correlation and many-electron 
effects in the system. 




Electron-ion collision energy (eV) 



FIG. 6: Comparison of the experimental Si IV DR spectrum 
(open circles) in the energy region — 6 eV with the final 
result of the MCDF calculation (solid line) . This energy range 
includes 2p e 3pn'l' (4 < n' < 6) and 2/ 3d 2 resonances. The 
nonresonant part of the recombination rate coefficient due to 
RR was subtracted from the measured curve (see text) 




Electron-ion collision energy (eV) 

FIG. 7: Comparison of experimental Si IV DR spectrum (open 
circles) in the energy region of the 2p 5 3s 31 n'V resonances 
with the results of our MCDF calculation (solid line). The 
nonresonant part of the recombination rate coefficient due to 
RR was subtracted from the measured curve (see text) 



In addition to calculations of the low-lying resonances 
associated with 3s — * 3p and 3s — * 3d core excita- 
tions, computations have been carried out also for all 
2p 5 3s 31 31' and 2p 5 3s 31 41' resonances which are found 
to occur starting at 69 eV. For this high-energy part of 
the spectrum, we expect only a rough agreement between 
our computations and experiment as the electronic struc- 
ture of the intermediate resonances now includes four 
open shells which do not allow any systematic enlarge- 
ment of the configuration basis. For this part of the spec- 
trum, therefore, the computations have been restricted 
to allow only one electron in the n = 4 shell. As seen in 
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Figure [7] only some of the resonances at 84 and 91 eV 
are reproduced by our simulations but they are shifted 
upwards in energy by about 1 eV. 

Apart from the radiative stabilization of the interme- 
diate resonances by El electric-dipole decay, we explored 
also the effects of higher multipoles (Ml, E2, M2) in the 
coupling of the radiation field. These 'multipoles' are 
well suppressed for neutral systems by at least 5 orders 
of magnitude but become important for highly-charged 
ions. For the initially triply ionized Si IV ions, these con- 
tributions are still negligible at the present level of accu- 
racy for the simulation of the DR spectra in Figures [5]-[7] 



coefficient for comparison with theoretical results (cf. 
Sec.lrVAl. 



When the electron-ion collision energy E is larger 
than the experimental energy spread AE, one can use 
q; dr/(2£ , /to ) 1 / 2 instead of ctdr for the convolution. 
When E < AE, the energy spread influences the out- 
come of the convolution of the cross section. In order to 
account for this effect, the low energy DR cross section 
was extracted by fitting 25 DR resonance line-shapes to 
the measured DR spectrum in the energy region below 
1.5 eV (Fig. [4}, independent of the theoretical predic- 
tions above (see Ref. |H] for a more detailed description 
of the method) . 



V. PLASMA RATE-COEFFICIENTS 

DR rate coefficients for a Maxwellian plasma can be 
derived from the experimental merged-beams recombi- 
nation rate coefficient and the theoretical cross section. 
This is done in the following sections. Of particular in- 
terest for astrophysical model calculations are the plasma 
DR rate coefficients in the electron temperature ranges 
where Si IV is formed in astrophysical plasmas. The ap- 
proximate temperature range where Si IV forms in pho- 
toionized and collisionally ionized plasmas can be ob- 
tained from the work of Kallman and Bautista [12] and 
Bryans et al. [43], respectively. For photoionized plas- 
mas, Kallman and Bautista [J2] find that the fractional 
Si IV abundance peaks at a temperature of 1 eV. The 
'photoionized zone' may be defined as the temperature 
range where the fractional abundance of a given ion ex- 
ceeds 10% of its peak value. For Si IV this corresponds 
to a temperature range of 0.8-1.4 eV. Using the same 
criterion and the fractional abundances of Bryans et al. 
|43j , for coronal equilibrium the Si IV 'collisionally ionized 
zone' is estimated to extend over a temperature range of 
4-10 eV. It should be kept in mind that these temper- 
ature ranges are only indicative. They depend, in part, 
on the accuracy of the underlying atomic data base. 



A. Derivation of the plasma DR rate coefficients 

The DR rate coefficient in a Maxwellian plasma is de- 
rived by convolving the DR cross section odr, with an 
isotropic Maxwell-Boltzmann electron energy distribu- 
tion as detailed by Schippers et al. [SJ |Hj. To derive 
a meaningful plasma DR rate coefficient from a total 
merged-beams rate coefficient there are some issues that 
require special consideration. 

Interference between DR and RR is typically unim- 
portant |45j . Hence here we subtract, the non-resonant 
RR contribution from the measured merged-beams re- 
combination rate coefficient. The applied merged-beams 
RR rate coefficient was the same as used for the cor- 
rection of the recombination rate coefficient [Eq. Q] at 
the reference point as well as that which was subtracted 
from the experimental merged-beams recombination rate 
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FIG. 8: Comparison between the experimental merged-beams 
DR rate coefficient (open circles with solid line) and the AU- 
TOSTRUCTURE calculation. The AUTOSTRUCTURE cal- 
culation was multiplied by a factor of 1.13 (see text). The 
dashed line is the theoretical result with account for the ex- 
perimental field ionization of high-n Rydberg states using the 
field ionization model of [Sj. The shaded area highlights the 
unmeasured purely calculated part of the composite DR rate 
coefficient. 

Field ionization of the loosely bound high Rydberg 
electron in the recombined ions can result from the mo- 
tional electric fields that the ions experience inside the 
storage ring bending magnets [BJ. The ion beam on its 
way to the detector passes the strongest electrical field 
in the dipole bending magnet in front of the detector. 
From this electric field the critical quantum number for 
field ionization is n cr i t = 17, i.e., in the present experi- 
ment only RR and DR involving capture into Rydberg 
levels with quantum numbers less than 17 contribute fully 
to the measured merged-beams recombination rate coef- 
ficient. Due to radiative decay of higher Rydberg states 
on the way from the interaction section to the detec- 
tor, the field ionization cutoff is not sharp but somehow 
smeared out ton > 17. Similar to the approach of Schip- 
pers et al. [HE] the missing DR resonance strength up to 
"max = 1000 was estimated from a theoretical calculation 



using the AUTOSTRUCTURE code [2S] . 



1000 



is an arbitrary upper limit beyond which no significant 
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contribution to the total DR cross section is expected. 
Although the AUTOSTRUCTURE code does not repro- 
duce the resonance structure below 5 eV (cf. Fig. 1 of 
[H]) as accurately as our present MCDF calculations, it 
reproduces the more regular structures of Rydberg res- 
onances between 5 eV and 8.4 eV when the calculated 
rate coefficient is multiplied by a constant factor of 1.13. 
The unmeasured DR contribution due to n > 17 exceeds 
the measured contribution by more than one order of 
magnitude. This is shown in Figure |U The DR contri- 
bution for AN = 3s -> 3p DR from n = 17-1000 was 
added to the measured spectrum by adding the differ- 
ence between the measured rate coefficient and the ad- 
justed AUTOSTRUCTURE result in the energy range 
8.39 eV — 9.01 eV (grey shaded area in Figure [8]). 

The 3s — ► 3d AN = series with its limit at about 
20 eV was not corrected for field ionization losses of 
high Rydberg states because its contribution is negligi- 
ble compared to the 3s — > 3p AN — scries limit. As 
the strengths of resonances contributed by AN > 1 DR 
fall much faster with increasing n than the strength of 
AiV = DR resonances, and because of the much smaller 
contribution of AN > 1 DR, the field ionized contribu- 
tion for AN > 1 DR with n > 17 was also not corrected 
for field ionization losses, either. 

The non-measured contribution to the plasma DR rate 
coefficient ranges in the photoionized zone from only 1% 
at 0.8 eV to 36% at 1.4 eV. It has its maximum (83%) at 
a plasma electron temperature of k^T c — 7.9 eV. In the 
collisionally ionized zone the contribution ranges from 
79% at 4 eV to 83% at 10 eV. The contribution falls 
off slightly towards higher temperatures and is still 72% 
at 1000 eV. The resulting plasma DR rate coefficient 
with and without the correction of field ionization losses 
is shown in Figure [9] 

For convenient use in astrophysical modeling codes the 
Si IV to Si III plasma DR rate coefficient atp 
using 



(DR.) 



was fitted 



< R) (T C ) = (T c )- 3 / 2 ^>exp(-£; 4 /fcBT c ). (9) 



i=i 



The fitting parameters Cj and are given in Ta- 
ble III The fit deviates by less than 1% from the 



experimentally-derived result in the temperature range 
0.01 eV - 10000 eV. 



B. Comparison with present theory 

In Figure [9] we compare the plasma DR rate coefficient 
derived from the MCDF calculation in the electron-ion 
collision energy range 0-6 eV and a plasma DR rate coef- 
ficient derived from DR merged-beams resonances mea- 
sured in the same electron-ion collision energy range. 
We find good agreement in the comparison of these two 
plasma DR rate coefficients. The plasma DR rate coeffi- 
cient generated by DR resonances, calculated by means 



TABLE III: Parameters for the fit of Equation (9| to the 
experimental plasma DR rate coefficient. Numbers in square 
brackets denote powers of 10. 



i 


c, (cm^- 1 ^ 2 ) 


Ei (eV) 


1 


2.13[-8] 


1.02[-2] 


2 


6.12[-8] 


5.00[-2] 


3 


1.10[-6] 


1.24[-1] 


4 


3.65[-6] 


2.44[-l] 


5 


1.45[-5] 


6.90[-l] 


6 


1.78[-5] 


1.53[+0] 


7 


3.05[-4] 


5.43[+0] 


8 


9.50[-3] 


8.81[+0] 


9 


1.89[-3] 


8.05[+l] 



of the MCDF method for electron-ion collision energies 
below 6 eV is somewhat lower than the rate coefficient 
generated from the experimentally-derived resonances in 
the same energy range. Between k-QT c = 0.01 eV and 
fceTc = 0.04 eV the agreement is better than 12%. Above 
a plasma electron temperature of fee^c = 0.04 eV the 
agreement is even better than 7%. 



C. Comparison with previous results 

In Figure [9] we also compare our experimentally- 
derived DR rate coefficient with the experimental result 
of Orban et al. [5] and with recent theoretical results of 
Gu (46] and Altun et al. [37]. 

In the collisionally ionized zone the rate coefficient of 
Gu |46j overestimates our experimentally-derived plasma 
DR rate coefficient by about 10%. In contrast, the 
plasma DR rate coefficient of Altun et al. 37J underesti- 
mates the experimentally-derived plasma DR rate coef- 
ficient in the collisionally ionized zone by about 10%. 
These deviations are within the experimental uncer- 
tainty. The difference between the two codes of about 
20% also gives an idea as to the uncertainty of these 
state-of-the-art DR calculations. 

In the photoionized zone the experimentally-derived 
plasma DR rate coefficient is decisively determined (be- 
tween 93% at 0.8 eV and 41% at 1.4 eV) by the 3s+e" -> 
3p Al' resonances which occur at electron-ion collision en- 
ergies below 1.5 eV. The calculations of Gu [IB] and 
Altun et al. [37] comprise DR for all AiV = chan- 
nels including the 3p4Z' resonances. Despite the convo- 
lution of the DR cross section with the plasma electron 
energy distribution one can recognize that the calcula- 
tions of Gu [IB] and Altun et al. [37] either underesti- 
mate the strengths of the 3p4/' resonances at low en- 
ergies or overestimate their energy positions. We can 
reproduce the qualitative trend of their plasma DR rate 
coefficients by shifting our experimentally derived 3p4Z' 
resonances to higher energies by about 0.15 eV. Both the- 
oretical plasma DR rate coefficients underestimate the 
experimentally-derived rate coefficient in the photoion- 
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FIG. 9: Experimentally derived Si IV to Si III DR rate coef- 
ficient in a plasma (thick solid line) comprising AJV = DR 
(Figs, [I] and [S), A TV = 1 and 2 DR (Fig. pj), and the the- 
oretical estimate for the unmeasured contributions of states 
with n > 17 for A TV = DR (Fig. [8}. The error bars denote 
the ±18% (la) experimental uncertainty in the absolute rate 
coefficient. The experimental results without DR extrapola- 
tion is shown by the thin solid line. Also shown are recent 
theoretical calculations of the DR rate coefficient by Gu |46] 
(thick dash-dot-dotted line, labeled Gu) and Altun et al. [37] 
(thick dash-dotted line, labeled Altun), and the recent exper- 
imental DR rate coefficient by Orban et al. [9] (thick dashed 
line, labeled Orban). The contribution from the experimen- 
tally measured DR resonances between and 6 eV is shown 
as the thin dotted line (labeled Exp < 6 eV) while the con- 
tribution of the MCDF calculation in the same energy range 
is shown as thick dotted line (labeled MCDF). A recent cal- 
culation of the plasma RR by Badnell [17] is shown as the 
thin dash-dot-dotted line labeled RR Badnell. The tempera- 
ture ranges where Si IV is expected to peak in abundance in 
photoionized plasmas (labeled PP) and collisionally ionized 
plasmas (labeled CP) are highlighted. 

ized zone. The rate coefficient of Altun et al. 37J is 
a factor of 0.71-0.86 lower than the experimentally de- 
rived rate coefficient. For the rate coefficient of Gu [IB] 
the factor is 0.83-1.02. 

The comparison of our experimentally derived plasma 
DR rate coefficient with the experimental result of Orban 
et al. |9J shows, that both data sets agree to within 25% 
in the temperature rage 0.1-1000 eV. In the photoion- 
ized zone the plasma DR rate coefficient of Orban et al. 
[9] is 18%-24% larger than our result. This larger de- 
viation than in the collisionally ionized zone is probably 
associated with differences in the data reduction process. 
In the collisionally ionized zone the plasma DR rate co- 
efficient of Orban et al. [9] is less than 5% lower than 
our result. The increasing deviation between the two 



plasma DR rate coefficients above ss 30 eV is probably 
attributed to the fact that, in contrast to the work of 
Orban et al. jS], the present experimentally derived DR 
rate coefficient also includes ATV = 1 and even ATV > 1 
DR. 

VI. CONCLUSIONS 

Electron-ion recombination of Si IV forming Si III was 
studied both experimentally using the merged-beams 
method at a heavy ion storage ring and theoretically by 
employing the MCDF method. We see good agreement 
in DR resonance strength and positions between the ex- 
periment and the MCDF calculations for ATV = DR 
in the investigated electron-ion collision energy range 0- 
6 eV. Below an energy of 1.4 eV the accuracy in the 
energy positions was better than 70 meV, the position of 
the resonance at the lowest energy was even accurate to 
within 20 meV. A great advantage of the MCDF method 
is that it is conceptually much simpler to implement than 
other many-body techniques and, hence, can be applied 
also to more complex shell structures — if enough com- 
putational resources are available. We currently plan to 
extend the code in order to make computations feasible 
for atoms and ions with (initially) two or three electrons 
in their valence shell. 

The present experimentally derived Si IV plasma DR 
rate coefficient agrees with the experimental result of Or- 
ban et al. [j]], to within the combined experimental errors. 
We found good agreement between the theoretical results 
of Gu [IB] and Altun et al. [37] and our experimental re- 
sult in the temperature range where Si IV forms in col- 
lisionally ionized plasmas. The agreement is reasonable 
at temperatures where Si IV is predicted to form in pho- 
toionized as. At temperatures below this the agreement 
becomes significantly worse with decreasing temperature. 
These findings demonstrate the necessity of benchmark- 
ing theoretical results with experiment, because modern 
theory still has difficulty calculating resonance energies 
reliably when the electron-ion collision energy and Ryd- 
berg level of the recombined ion are small. 
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